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Abstract 

We report the discovery of the afterglow of the X-ray rich, long-duration gamma-ray 
bm'st GRB 011211 and present evidence for oscillatory behaviour in its early optical 
light curve. The time-scale of the fluctuations, ~1 hour, is much smaller than the 
time of the observations, ~12 hours from the onset of the gamma-ray burst. The 
character and strength of the fluctuations are unprecedented and are inconsistent 
with causally connected variations in the emission of a symmetric, relativistic blast 
wave, i.e. flux variations which are produced uniformly throughout the shell sur- 
face are ruled out. Therefore, the wiggles are the result of spherically asymmetric 
density or energy variations. Additionally, there is evidence for fluctuations in the 
X-ray afterglow light curve. If real, the resulting difference in the observed time 
of the peaks of the short-term variations at X-ray and optical frequencies, would 
demonstrate that the energy content across the jet-emitting surface is not uniform. 
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1 Introduction 



GRB 011211 was detected on 2001 December 11.798 UT (|Gandolf]i l2001al ) 
with one of the Wide Field Cameras (WFCs) on board the Italian-Dutch 
satellite BeppoSAX. The prompt gamma-ray emission lasted ~270 s, making it 
one of the longest bursts o bserved with the s atellite. Followi ng the distribution 
of a 2' radius error circle ( Gandolfi . 2001b[ ). we identified ( Gray et al. . 2001) 
the optical counterpart of the gamma-ray burst (G RB). An X-ray afterglow 
subsequently detected (jSantos-Lleo et al. [ l200lh and its redshift 



was 1 



measurea 



yia absorption lines in the optical spectrum to be ^ = 2.140 (IFruchter et al 



2M 



Holland et al. , 2002[ ) . The analysis of the X-ray spectrum (jR.eeyes et al 



20021 ) showed emission lines arising in metal-enriched material with an outflow 
yelocity of ~0.1c. The presence of such line features strongly sug gests a massiye 
stellar progenitor, but the details remain model dependent (jLazzati et al 
2002ah . 



Afterglow obseryations of long-duration GRBs typically show their flux to 
decline as a power-law in time. Howeyer, recent continuous early-time mon- 
itoring of GRB optical light curyes h as shown eyidence for yariations or 
bumps superp o sed o n power-law decays ( Bersier et al. . 20031 : Fox et al. . 20031 : 
Uemura et all l2003h . We note that these are not supernoya bumps as ob- 
seryed in the optical afterglow (OA) light curyes of sev e ral bu rsts (see e.g. 



Bloom et al.1 Il999l : iGarnayich et al.L 120031 : iHiorth et all 12003^ We present 



here for the first time a comparative study of multi- wavelength oscillations ob- 
served in the light cu rve of GRB 1 1211. Despite intensive monitoring of many 
GRB counterparts (IHiorth et al 1 I1999I: IStanek et al.l [ I1999I: iHalpern et al 



Gorosabel et al. 



2OO3I ). 



2OO3I : iBurenin et all 120031 : iLaursen fc Stanekl . l2003l : 

the only other burst for which short-term time-scale (less than one hour) vari- 
ations have bee n detected in its opti cal light curve is GRB 021004, albeit in a 
random fashion ( Bersier et al. . 2003| ). The oscillatory behaviour and the rapid, 
sharp decline and rise of the pulses of the light curve of GRB 011211 presented 
here, are unlike any other light curve reported before. Our observations pro- 
vide a useful probe of the GRB explosion, the structure of the emerging jet, 
as well as the structure of the surrounding environment. 
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Fig. 1. Images of the field of the optical transient of GRB 011211. Left: The R-band 
DSS-2 image. Right: The discovery R-band image of the GRB 011211 optical after- 
glow (OA), taken with the N OT. The position of the OA is marked with an arrow. 
The BeppoSAX localization \Ga,ndoln 
Both images cover a square of 5' x 5'. 



2001 A ) is indicated by a circle of radius 2'. 



2 Observations &; data analysis 



We observed GRB 011211 ~9.6 hours after the GRB trigger (starting 2001 
December 12.22 UT) with the Nordic Optical Telescope (NOT) at the Canary 
Island of La Palma and with the Danis h 1.54-m telescope at La Silla in Chile, 
in standard B-, V-, R-, and /-bands f Grav et al. . 200l[ ). A detailed log of 



our observations is listed in Table 1. A comparison of our first image with 
the Digitized Sky Survey 2 (DSS-2) frames, identified a new source shown in 
Fig. 1. Subsequent observations with the NOT revealed that the source was 
fading, thus establishing that the nature of the transient was consistent with 
being the OA of GRB 011211. 

Reference stars for astrometry were collected using the large field of view 
provided by the Danish Faint Object Spectrograph and Camera (DFOSC) 
mounted on the Danish 1.54-m (13.7 x 13.7 arcmin^). Based on ~50 USNO 
stars per image, we determined 11 independent positions of the afterglow. The 
mean value of the 11 afterglow coordinates is RA(J2000) = ll'^15™'17f98 and 
Dec(J2000) = -21°56'56'.'l with an error of 0'.'3. 

We obtained very well sampled early i?-band observations of the OA ~10-14 
hours after the burst. Our data indicated significant tempo ral variations su- 



perpo sed on the power-law fit to the light curve of the GRB (|Jakobsson et al 



2OO3I hereafter J03). To enhance the visibility of these features, we subtracted 



from our data the power-law decay contribution (ai = —0.95) in the optical 
flux; the result is shown in Fig. 2 (top panel). The deviation of the data points 



3 



Table 1 

Photometry of the afterglow of GRB 011211. The magnitude of the optical afterglow 
was measured relative to eight stars in the field. The calibrated magnitudes of these 
stars are given in JOS. Due to saturation we only used a subset of these stars in 
some of our images. The magnitudes were calculated using point-spread function 
fitting photometry. 



Date (UT) 


Tel. 


Magnitude 


Seeing 


Exp. time 


(2001 Dec) 






(arcsec) 


(s) 


B-band: 










12.2816 


NOT 


21.229 ± 0.047 


1.1 


300 


12.3539 


1.54-m 


21.262 ± 0.089 


1.4 


200 


V-band: 










12.2907 


NOT 


20.761 ± 0.080 


1.0 


300 


12.3586 


1.54-m 


20.987 ± 0.042 


1.2 


300 


R-band: 










12.2181 


NOT 


20.012 ± 0.033 


1.7 


300 


12.2227 


NOT 


19.877 ± 0.029 


1.4 


300 


12.2272 


NOT 


19.918 ± 0.022 


1.3 


300 


12.2831 


1.54-m 


20.368 ± 0.040 


1.5 


600 


12.2912 


1.54-m 


20.298 ± 0.033 


1.5 


600 


12.2994 


1.54-m 


20.230 ± 0.038 


1.4 


600 


12.3075 


1.54-m 


20.285 ± 0.042 


1.4 


600 


12.3157 


1.54-m 


20.274 ±0.051 


1.3 


600 


12.3238 


1.54-m 


20.330 ± 0.046 


1.1 


600 


12.3633 


1.54-m 


20.539 ± 0.063 


1.1 


300 


12.3681 


1.54-m 


20.616 ± 0.100 


1.0 


300 


I-band: 










12.2861 


NOT 


19.965 ±0.110 


0.9 


300 


12.3350 


1.54-m 


19.951 ± 0.054 


1.3 


300 


12.3397 


1.54-m 


19.981 ±0.071 


1.8 


300 


12.3491 


1.54-m 


20.064 ± 0.061 


1.4 


300 
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Fig. 2. Comparison of the residuals in the optical and X-ray light curves of the af- 
terglow of GRB 011211. Top : The R-band data are drawn from Table 1 and the lit- 
erature \Holla.nd et a,l\ . \200^ }. Middle: The same as in the top panel, except here we 
also plot the I-, V-, and B-band residuals (data taken from Table 1). Offsets relative 
to the R-band were applied using color information from the spectral energy distri- 
bution (JOS). Bottom: The XMM-Newton X-ray light curve residuals (0.5-10 keV). 
In order to compare the X-ray results with the optical, the units are expressed in 
terms of —2.5 x log{fiux). In all three panels, the horizontal error bars are equal to 
the integration time for each data point, and t is the time since the burst. 
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from a pure power-law decline is significant (>5(j) as derived from Spearman's 
rank correlation. The residual weighted mean is 0.01 ± 0.03 with a root mean 
square (rms) of 13% (after correcting for the photon noise). The amplitude 
of the largest oscillation is ~0.25 magnitudes. We detect two maxima with a 
single minimum in-between and by fitting Gaussians to the data set we mea- 
sure the time-scale of the temporal oscillations to be At^ = 0.9 ±0.2 hours. In 
the middle panel of Fig. 2 we display the combined multicolor light curve of 
GRB 011211; a wave- like behaviour similar to that in the i?-band light curve 
is evident in the combined /-, V-, and i?-band residuals. 

The bottom panel of Fig. 2 shows the temporal variations present in the resid- 
uals of the X-ray afterglow of GRB 0112 11 observed with t he X- ray Multi 
Mirror- Newton satellite (XMM- Newton) (|Santos-Lleo et all I2OOIL The re- 



duction of the X-ray data set is described in detail in J03. We observe sig- 
nificant deviations from the power-law decline (ax = —1.62) in the X-ray 
flux (>3(T as derived from Spearman's rank correlation). Although the X-ray 
fluctuations are not as signiflcant as the optical ones, the value of the (X- 
ray) Spearman rank correlation coefficient translates to a probability of less 
than 3 x 10"'^ that this correlation is a statistical fluctuation. The weighted 
mean of the residuals is 0.00 ± 0.19, with an ~11% rms variation around 
the mean flux (after correcting for the photon noise). There is weak evidence 
that the fluctuations are of a similar oscillatory nature as the optical ones. 
The fluctuation time-scale is comparable to that observed in the optical, with 
Atx/Atfi = 0.75 ±0.20. From the short overlap of the optical and X-ray data, 
we flnd that the oscillations are not in phase with a time difference of at least 
0.70 ± 0.15 hours (as derived from the cross correlation function). 



3 Discussion 



The short-term wave-like behaviour of the GRB 011211 optical light curve 
reported here is unprecedented. This behaviour was not detected in th e earlier 
analysis of the OA of GRB 011211 reported bv lHoUand et al.l (|2002h . due to 
insufficient data sampling. Moreover, the vari ations detected in the optical 
afterglow of GRB 021004 ( Bersier et al. . 2003| ) appeared highly random and 
erratic, and all other GRBs have displayed variations with much longer time- 
scales. We discuss below the origin of the bumps in the optical light curve of 
GRB 011211. 



The simplest afterglow model, where a relativistic jet decelerates as it ex- 
pands into the ambient matter leading to a radiative output with a char- 
acteristic power-law decay, has been remarkably succ essful so far in accom- 
modating the present data ( van Paradiis et al. . 2000[ ). How does this model 
account for the unusual optical light curve of GRB 011211? The strong tem- 
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poral variations in the early afterglow could be interpreted as a result of i) re- 
freshed shocks created as the leadi ng edge of the iet decelerates and is caught 
up by slower-moving iet gas (e.g . Rees fc Meszarosl . 19981 : Panaitescu et all 
19981 : Ramirez- Ruiz et al. . 2001lJ ). ii) the relativistic jet impacting an exter- 



nal medium of variable dens i ty (e.g. Wang fc Loebl. 2000l: Rarnirez-Ruiz et al. 



2nnial : iLazzati et al.l . liocml iHevl fc Pernal. 120031: iNakar et al.l. 120031) . or mj a 



non-uniform jet structure (e.g. iMeszaros et al. . Il998l : iKumar fc Piran . 2000| ) 



We note that a smooth jet traveling through a clumpy medium would quickly 
cease to be homogeneous, i.e. it is possible that ii) could give rise to Hi). 
Refreshed shocks covering a large fraction of the emitting surface can only 
increase the energy of the blast wave, and therefore cannot account for the 
rapid decay seen in the optical data (top panel of Fig. 2). 



The sharpness of the features observed at optical frequencies is difficult to 
reconcil e with a density disco ntinuity covering most of the visible surface of 
the jet ( Nakar fc Piran . 2003[ ). On the other hand, the short time-scale of the 
oscillations provides interesting upper limits of ~1-10 AU on the size of the 
clumps around the source. These limits are lower than the flu ctuation ampli- 
tudes seen on similar scales in the local interstellar medium ( Wang &: Loebl 
2000l ). though they may reflect the length scale of comet-like clumps obser ved 



in ring nebulae surrounding massive stars (e.g. Garcia-Segura et al. . 1996| ). 



A promising alternative interpretation is energy variations within the expand- 
ing jet. It is possible for a large fraction of the emitting surface of the jet 
to become active and the flare region to remain small, because the energy 
content of the jet varies strongly as a function of angle. At the time t of the 
observations the region contributing to the total observed flux would have a 
transverse size of 2rct, where F is the Lorentz factor of the jet. Only regions 
the size of FcAt can produce afterglow fluctuations with time structures of 
At ^ t. In other words, the jet-emitting surface is peppered with many regions 
(hot spots) on a scale much smaller than the narrow cone visible along the line 
of motion (see Fig. 3). Variations in the initial conditions as a function of the 
opening angle could spread the causally disconnected regions out along the fine 
of sight to the observer such that the emission would arrive at different times. 
The corresponding overall emission at a given wavelength is then averaged 
over the observed region, and it may increase or decrease locally depending on 
the number, intensity, and spectrum of the different flare sites across the jet 
surface. This can happen either because the radiating material is in narrow 
flngers or, as the jet slows down, we observe afterglow emission not just from 
material moving directly towards us, but also from a wider range of angles. 
Thus, an inhomogeneous energy jet is a viable explanation for the observed 
rapid variations in the light curve, with the size of a hot spot determining the 
duration of the bump. 

From our optical observations alone we can safely conclude that the wiggles 
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Fig. 3. GRB tomography. The emission from a radiating object moving close to the 
speed of light is beamed into a very narrow cone, (a) A relativistic jet that radiates 
uniformly would have an emitting region with a transverse size of2Tct. The afterglow 
profile would be the sum of many emitting regions, and thus, should be smooth. 
(b) Only regions the size of FcAt can produce afterglow fluctuations with a time 
structure of At. The same jet pattern observed face-on is shown in (c). The photon 
flux emitted per unit solid angle and per unit time at this flare site should be larger 
than that arising from the rest of the blast wave. The corresponding emission would 
then average over the observed region and it may increase (or decrease) depending 
on the intrinsic properties of different flare sites across the jet surface. For example, 
a single region with a T that is ten times larger than for other regions could produce 
emission that doubles the overall count rate in a different energy band. 
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in the light curve of GRB 011211 are the result of spherically asymmetric 
density or energy variations, i.e. variations that cover less than the observed 
1/r region. The question remains: Can we decide between the above two 
scenarios. In principle, simultaneous multi-wavelength observations can help 
settle this question (see Salmonson & Ramirez-Ruiz, in preparation, for a de- 
tailed treatment) and, although the X-ray variations are of lower significance 
than the optical, in what follows we explore the implications of fluctuations in 
both bands. As shown by JOS, the X-ray band is located above the cooling fre- 
quency at the time of our observations, where density variations of the ambient 
medium hav e a mu ch weaker impact on the emitted flux (see e.g. equation 7 in 
Nakar et al. . 2003| ). Given the observed ~11% variation in the X-ray flux, we 



would then expect a much stronger correlated optical one, which we do not ob- 
serve. On the other hand, uncorrelated fluctuations at different energy bands 
can be produced by hot spots (small regions with higher energy and/or F). 
Since F strongly affects the synchrotron peak frequency, a situation can easily 
arise where we observe a hot spot releasing most of its energy in the X-rays, 
without a significant optical contribution above the continuum. Conversely, 
an optical/infrared bump can be detected without an X-ray counterpart, at 
least when the cooling frequency of the afterglow continuum spectrum is lo- 
cated between the optical and the X rays. Thus, the combined departure from 
a power-law decay observed in the optical and X-ray light curves seems to 
disfavour a scenario in which the density of the ambient medium varies and 
the energy is kept constant across the jet surface. 

The above reasoning assumes that the X-ray and optical photons are emit- 
ted from the surface of the expanding jet, an assumption justified by the fact 
that the variations seen at different energies have similar time-scales. Inter- 
estingly, the analysis of the X-ray spectrum during this time frame uncovered 
the presence o f emi ssion line features arising from metal-enriched material 
( Reeves et al. . 2002| ). It is feasible that a fraction of X-ray photons reach 



the observer after be ing reflected by material away from the line of sight 
( Lazzati et al 1 l2002at) . While the presence of such spectral features is not 



conclusive (but see e.g. Watson et al. . 2003| ). the X-ray wiggles observed here 



could be attributed to variations in the line emission intensity, which in turn 
are caused by changes in the illuminating continuum. The X-ray afterglow 
could then, at least in part, be due to the continuing power output from the 
GRB central engine inter acting with the envelope of a massive progenitor star 
( Rees fc Meszarosl . l2000l) . 



CoUapsar ( MacFadven et al. . 2001 ) or mag netar-like (|Wheeler et all I2OOOI ) 



models not only provide a natural scenario for a sudden burst followed by a 



^ This is because the observed time-scale of the variations (~1 hour) is much shorter 
than the overall elapsed time after the burst trigger (^^12 hours), when angular 
smoothing would have smeared out the light curve variations. 
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decaying energy output but could also easily imprint a non-unif orm structure 



to th e GRB jet as it bores its way through the stellar mantle (jZhang et al. 



2003| ). Further data on early X-ray spectral and temporal features provided 



by the Swift satellite will offer additional clues to the nature of the progenitor 
and the relativistic jet that it triggers. 
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